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Abstract 


The  scattered  fields  from  a  spherical  body  eccentrically  located  within  a  host 
sphere  are  found  by  a  derivation  that  satisfies  the  boundary  conditions  at 
both  interfaces.  The  source,  which  may  be  composed  of  any  linear  combina¬ 
tion  of  S  and  P  waves,  is  also  arbitrarily  located  within  the  host  sphere.  The 
scattering  system  has  applications  in  seismic  scattering,  since  the  scattering 
interaction  between  a  scatterer  and  the  Earth’s  surface  is  significant  when  the 
scatterer  is  located  near  the  surface.  This  interaction  can  affect  subsequent 
seismograms. 
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Figures 


1  Geometry  of  scattering  system  showing  coordinate  systems  of 
host  (xi,  yi,  2i),  scatterer  {x2,  y2,  -22),  and  source  {xz.yz,  Z3)  .  . 


1.  Introduction 


Scattering  of  elastic  waves  by  spherical  bodies  is  often  studied  in  the  con¬ 
text  of  seismic  theory.  The  simplest  system  for  which  an  analytical  solution 
can  readily  be  obtained  for  the  entire  scattered  elastic  wave  field  is  that  of 
a  sphere  in  an  infinite  elastic  medium  (Ying  and  Truell,  1956;  Einspruch, 
Witterholt,  and  Truell,  1960;  Dubrovskiy  and  Morochnik,  1989).  More  com¬ 
plicated  techniques  must  be  employed  to  take  into  consideration  the  effects 
of  multiple  scattering  (Wu,  1985;  Wu  and  Aki,  1985)  and  irregularities  in 
the  geometry  of  scatterers  (Waterman,  1969,  1976;  Varatharajulu  and  Pao, 
1976).  When  an  inhomogeneity  is  located  near  an  interface,  the  interaction 
between  the  inhomogeneity  and  the  interface  is  not  insignificant  compared 
to  the  incident  field,  and  must  be  considered  when  the  scattered  fields  are 
calculated. 

In  this  report,  we  derive  the  equations  describing  the  scattered  field  from 
a  spherical  scattering  body  eccentrically  located  within  a  host  sphere.  The 
interface,  in  this  case,  is  the  host  sphere  (free)  surface,  which  can  be  made 
large  compared  with  the  size  of  the  scatterer.  Such  a  system  is  of  interest 
in  seismology,  since  it  can  be  used  to  describe  the  scattered  fields  due  to 
inhomogeneities  located  near  the  surface  of  the  Earth  (host  sphere).  In  order 
to  solve  this  problem,  we  expand  the  relevant  wave  fields  in  terms  of  vector 
spherical  harmonics  and  use  the  interface  boundary  conditions  at  both  in¬ 
terfaces.  In  order  to  solve  both  sets  of  boundary  conditions  simultaneously, 
we  translate  the  spherical  vector  harmonics  representing  the  individual  fields 
between  the  coordinate  systems  of  the  host  and  scattering  spheres.  This  work 
is  similar  to  previous  works  in  which  eccentric  sources  are  embedded  within 
a  spherical  host  (Thompson,  1973;  Glenn  et  al,  1985;  Rial  and  Moran,  1986; 
Zhao  and  Harkrider,  1992);  however,  in  addition  to  an  eccentric  source,  we 
also  include  an  eccentrically  located  scatterer.  The  notation  used  in  this  re¬ 
port  is  similar  to  that  of  a  previous  paper  on  light  scattering  firom  a  host 
sphere  containing  an  eccentric  inclusion  (Videen  et  al,  1995). 
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2.  Theory 


The  geometry  of  the  scattering  system  is  shown  in  figure  1.  A  host  sphere 
of  radius  a  and  elastic  moduli  A  and  fx  is  centered  on  the  Xi,  yi,  Zi  coordi¬ 
nate  system.  A  scattering  sphere  of  radius  a'  and  elastic  moduli  A'  and  n' 
is  centered  on  the  X2,  y2)  -^2  coordinate  system,  located  at  (xi  =  0,  yi  —  0, 
zi  =  d).  A  source  is  centered  on  the  Xz.ys,  Zz  coordinate  system,  located  at 
(n  =  Tj,  61  =  Os,  <^i  =  0)  such  that  Zz  is  parallel  to  fg,  and  yi  is  parallel 
to  yz-  A  time  dependence  of  exp(— iwt)  is  implicit.  Real  seismic  sources  have 
a  more  complicated  time  dependence.  For  linear,  homogeneous  systems,  we 
can  solve  the  equations  for  a  general  time-dependent  source  using  integral 
Fourier  transformations: 


= 

poo 

/  F{kj)exTp{—iXjkj)dkj, 

J  —  00 

(1) 

1 

(2) 

F(kj)  = 

In  this  report  we  solve  for  the  amplitudes  of  the  scattered  fields  at  specific 
frequencies  F  (kj),  which  we  can  transform  to  find  the  time-dependent  am¬ 
plitudes  using  equation  (1). 


Figure  1.  Geometry  of  scattering  system  showing  coordinate  systems  of  host  (xi,j/i,2i), 
scatterer  (x2,j/2,^2),  and  source  {xz,yz,zz). 
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We  reach  the  general  solution  by  satisfying  the  boundary  conditions  on  both 
spherical  interfaces  simultaneously.  In  doing  so,  we  consider  the  fields  inci¬ 
dent  on  the  interfaces  of  the  two  subsystems  (host  and  scattering  spheres) 
separately.  We  use  the  vector  spherical  harmonics  (Hansen  vectors),  which 
have  the  following  form  in  this  derivation: 
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where  the  prime  denotes  a  derivative  with  respect  to  the  argument,  the  index 
j  corresponds  to  the  coordinate  system  used  {j  =  1,2),  and  z^\krj)  are  the 
spherical  Bessel  functions  of  the  first,  second,  third,  or  fourth  kind  (p  = 
1,2, 3, 4)  (Ben-Menahem  and  Singh,  1981),  and 


P™(COS0,)  = 


(2n  -f  l)(n  —  m)\ 
■M  2(n-t-m)! 


Pr(cos0,)  , 


(6) 


where  P^{cos6j)  are  the  associated  Legendre  polynomials. 

It  is  convenient  to  consider  the  vector  surface  harmonics,  Pnm?  arid 
Cnm  when  solving  the  boundary  conditions.  The  vector  surface  harmonics 
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are  defined  in  terms  of  the  vector  spherical  harmonics: 
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We  start  the  derivation  by  expanding  the  fields  in  terms  of  the  vector  spheri¬ 
cal  harmonics  and  satisfying  the  boundary  conditions  on  the  interface  of  the 
scattering  sphere. 


2.1  Scattering  Sphere 

First  we  examine  the  fields  that  strike  the  outer  surface  of  the  scattering 
sphere.  We  consider  an  arbitrary  field  incident  on  the  system  that  can  be 
expanded  with  the  spherical  Bessel  functions  of  the  first  kind,  jn{kr2)  : 


u: 
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where  k^  and  kp  are  the  compression  and  shear  wavenumbers  external  to 
the  scattering  sphere,  respectively.  Similarly,  k[^  and  k'^  are  the  compression 
and  shear  wavenumbers  internal  to  the  scattering  sphere,  respectively.  The 
coefficients  44’  44  represent  the  field  from  the  source.  Similarly, 

the  scattered  field  from  the  scattering  sphere  may  be  expanded  with  the 
spherical  Bessel  functions  of  the  third  kind,  /i^^)(/cr2): 
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This  scattered  field  interacts  with  the  host  sphere  (free)  surface,  which  reflects 
incoming  spherical  waves  onto  the  scattering  sphere.  These  may  be  expanded 
with  the  spherical  Bessel  functions  of  the  fourth  kind,  h^^kr2) : 

OO  n 

/  V  /  V  ^nm^nm,2  ^  ^nm^^nTn,2  •  ^mn^^nm,2 
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IV 
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The  fields  internal  to  the  scattering  sphere  may  be  expanded  into  standing 
waves  with  spherical  Bessel  functions  of  the  first  kind: 

OO  n 

“?n.  =  E  E  ^liL™ ,,  +  ,,  +  dS,NW ,,  (13) 

n=0  m=—n 
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With  the  expansions  of  the  fields  defined,  we  can  now  satisfy  the  boundary 
conditions  at  the  interface  of  the  scattering  sphere.  The  stress  vector  across 
a  spherical  surface  for  a  displacement  field  given  by 

u  =  Uir)Pnm  +  Vir)Bnm  +  W{r)Cnm  (14) 

can  be  calculated  as  (Ben-Menahem  and  Singh,  1981) 
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In  order  for  the  displacement  field  and  the  stress  tensor  to  be  continuous, 
the  following  six  boundary  conditions  must  be  satisfied  at  r2  =  a' : 

{kar2,  kpr2)  +  {kar2,  kp'^2)  +  Cnm  kpr2)  =  D!;^^  (k'a'<'2,  k'pr2)  ,  (16) 

^nm  {kar2,  kpr2)  +  {kar2,  kpr2)  +  {kar2,  kpr2)  =  [k'^r2,  k'pT'^  ,  (17) 

{kpr2)  +  {kpr2)  +  ( V2)  =  (/c^ra)  ,  (18) 
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We  are  interested  primarily  in  the  external  fields  of  the  scattering  sphere. 
By  substituting  equations  (16)  to  (18)  into  equations  (20)  to  (21),  we  can 
eliminate  the  internal  field  coefficients  and  express  the  scattered  external  field 
coefficients  in  terms  of  the  incoming  field  coefficients.  Then,  by  expanding  the 
vector  surface  harmonics  in  terms  of  the  vector  spherical  harmonics  (given 
in  eq  (10)),  we  can  derive  the  following  expressions  for  the  vector  spherical 
harmonic  coefficients: 

i  (22) 

nm  nm^nm  '  nm^nm  '  nm  nm  '  nm  nm  ’  v  / 

J,(2)  _  „(2)  q22  ,  (2)  d22 

and 

^(3)  =.(1)031  +a(3)Q33  I  c(l)/?31  +£(3)^33  (24) 

Explicit  expressions  for  Qi^,  Q3i^,,  Rll^, 

and  i?3^  are  given  in  appendix  A. 


2.2  Host  Sphere 

We  next  examine  the  fields  in  the  host  (Earth-centered)  coordinate  system. 
First  we  consider  a  source  that  can  be  expanded  with  spherical  Bessel  func¬ 
tions  of  the  third  kind,  hl^^{kri): 

oo  n 

„1  _  V  V  -<-p(3)]vj{3)  /OS') 

n=0  m=—n 
oo  n 

=  EE  i;W(n)P....,,,+f;2(ri)B„„,,  +  £i2(r,)C.„,,  . 
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The  scattered  and  incident  fields  from  the  scattering  sphere  may  be  expanded 
with  the  spherical  Bessel  functions  of  the  third  kind,  h^^\kri), 
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(27) 


and  the  spherical  Bessel  functions  of  the  fourth  kind,  hj^^(kri), 

oo  n  ,  . 

^ref  ■”  zL/  Z-/  ffnm^nm,!  ‘  9nm^^^nm,l  ^  9nm^^nm,l 
n=0  m=—n 
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In  appendix  B,  we  provide  expansions  whereby  relationships  can  be  derived 
between  the  and  and  and  and  coefficients.  The 

expansions  (26)  and  (27)  are  valid  in  the  region  where  n  >  d:  i.e.,  near  the 
host  sphere  (free)  surface.  Hence,  we  need  not  be  concerned  with  the  pole 
that  occurs  at  ri  =  0  for  the  scattered  and  incident  fields,  since  the  expansion 
is  not  valid  in  this  region.  Since  the  stress  vector  across  a  firee  surface  must 
be  zero,  the  following  boundary  conditions  apply  at  ri  =  a: 

(A  +  2,<)  +  -]  (ESLin)  +  F™(r,)  +  G<,«  (r,))  =  (28) 

ari  ^  J  ^ 

n{n  +  1)^  (F^m(7i)  +  Ej,^{ri)  +  G),^^{ri))  , 

a±  (£S;(n)  +  FS(n)  +  G®  (n))  +  (29) 

(E™  (ri)  +  E,SJ^(ri)  +  GSJ,(ri))  -  (E£(ri)  +  E®(ri)  +  GjS,(ri))  =  0  , 

^  (ES(rO  +  F®(r,)  +  G2(r.))  -  E®  (r.)  +  F^in)  +  G2(r,)  =  0  .  (30) 

As  with  the  scattering-sphere  coefficients,  we  express  the  coefficients  in  the 
host-sphere  coordinate  system  in  the  following  form; 

o(i)  -  (eW  I  fW\  5I1  4.  ,  (31) 

9nm  \^nm  '  J  nm )  ^nm  '  \^nm  ~  J  nm )  ^nm  ’  v  ) 

o(2)  =  fe(2)  .  f{2)\  522  (32) 

ynm  \^nm  ^  Jnmj  ‘^nm  )  \  > 

and 

o(3)  _  (M)  .  fW\  531  ,  (  (3)  ,  fi3)\  533  (33) 

ynm  ~  \^nm  ^  JnmJ  ^nm  “  \^nm  ~  Jnmj  ‘-'nm 

In  appendix  A,  we  give  explicit  expressions  for  and  5^^. 


2.3  Fields  Interior  to  the  Host  Sphere  (Earth) 

The  displacement  fields  given  by  equations  (10)  to  (12)  are  equivalent  to 
the  displacement  fields  given  by  equations  (25)  to  (27).  By  performing  a 
translation  of  the  spherical  vector  harmonics  from  the  host  coordinate  system 
to  the  scattering-sphere  coordinate  system  (using  the  results  of  app  B),  we 
can  put  equations  (31)  to  (33)  into  the  following  form: 
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Substituting  equations  (22)  to  (24)  into  these  expressions  yields 
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(38) 


fV)  = 


and 


OO 

E  (l)  ((^Z\  c22  r>(m,l)\  ,  (1)  ( p31  c-22  ry{m,l)\  , 

J  '  ^n'm  *^71771*^ n',n  J  ' 


The  and  coefficients  are  related  by  equations  (B-14)  to  (B-16) 
(app  B).  Since  the  coefficients  are  the  coefficients  for  the  incident  field, 
which  is  assumed  to  be  known,  the  can  be  easily  determined,  and  equa^ 
tions  (37)  to  (39)  represent  three  sets  of  equations  containing  three  sets  of 
unknowns  Numerical  results  based  on  this  analytic  solution  can  be  ac¬ 
quired  following  the  same  algorithms  used  in  light  scattering  (Videen  et  al, 
1995).  Equations  relating  the  various  fields  are  given  by  equations  (22)  to  (24) 
and  (31)  to  (33);  the  coefficients  used  in  these  equations  can  be  calculated 
from  the  equations  given  in  appendix  A.  Computational  results  obtained 
from  similar  expressions  in  the  field  of  light  scattering  have  shown  that  the 
coefficients  of  the  scattering  harmonic  terms  tend  to  have  a  negligible  con¬ 
tribution  to  the  scattered  fields  for  n  >  ka  +  4  (ka)^^^  -I-  2  (Wiscombe,  1980; 
Bohren  and  Huffman,  1983).  It  is  likely  that  a  similar  truncation  criterion 
holds  for  the  elastic  scatter.  With  the  truncation  of  the  summations  given 
by  equations  (37)  to  (39),  a  solution  can  be  found  through  matrix  inversion. 
With  the  and  known,  can  be  determined  by  equations  (22) 
to  (24).  Hence,  the  entire  displacement  field  within  the  host  sphere  can  be 
found. 
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2.4  Incident  Field 


Although  the  solution  derived  in  section  2.3  is  completely  general,  in  many 
applications  we  are  interested  in  a  specific  case  where  a  point  source  is  located 
at  some  position  within  the  host  sphere  (ri  =  =  0).  We  specify 

a  source  coordinate  system  {xz,  yz,  Zz)  such  that  Zz  is  parallel  to  r*  and  yi 
is  parallel  to  yz,  as  shown  in  figure  1.  The  derivation  is  valid  only  when  the 
source  is  closer  to  the  center  of  the  host  sphere  than  is  the  scattering-sphere 
center  {vg  <  d).  The  limitation  of  the  validity  is  a  result  of  the  expansion 
that  we  use  in  translating  the  vector  spherical  harmonics.  If  the  source  is 
farther  from  the  center  of  the  host  sphere  than  the  scattering  sphere  (r^  >  d) , 
then  we  would  need  to  use  a  different  set  of  translation  coefficients,  and  the 
derivation  would  be  shghtly  different  (Stratton,  1941).  The  source  can  be 
expanded  as 

oo  n 

=  E  E  '•SlS.s  +  +  aS.nS.,3  .  (40) 


This  field  can  be  expressed  in  a  coordinate  system  whose  origin  is  the  same  as 
the  host  (Earth)  coordinate  system,  but  whose  axes  are  parallel  to  the  source 
coordinate  system  by  a  translation  of  a  distance  Vg  in  the  —zz  direction.  The 
source  in  this  (x4,  y4, 24)  is 


OO  n 

n=0  m——n 

We  can  rotate  this  vector  field  using  the  results  of  appendix  B  to  give  the 
coefficients  of  equation  (25). 

We  provide  formulae  by  which  the  incident  field  coefficients  can  be  derived 
for  two  cases.  These  represent  simple,  but  solvable  sources.  The  incident  field 
coefficients  from  more  complicated  sources  can  be  substituted  for  the  coeffi¬ 
cients  given  in  this  section.  In  the  first  case,  we  consider  the  field  due  to  an 
explosion.  In  this  case,  the  source  pressure  wave  is  not  angularly  dependent, 
and  we  can  express  it  as  — or  as 


AiS,  =  AiSf  =  0  , 

a£  =  Aifi"  =  0  ■ 


(42) 

(43) 

(44) 


11 


In  the  second  case,  we  consider  a  source  that  emits  shear  waves  only,  whose 
shearing  stress  =  To  sin  63  is  applied  at  the  inner  surface  of  a  spherical 
cavity  centered  on  the  source.  In  this  case  the  source  coefficients  are  of  the 
form 


0 

II 

11 

(45) 

(46) 

^(3)  _  ^(3)s  _  Q 

(47) 

If  the  shear  is  in  another  direction,  we  can  use  the  vector  rotation  results 
given  in  appendix  B.  The  final  expressions  for  the  source  coefficients  in  the 
host  sphere  coordinate  system  are  given  by  equations  (B-31)  to  (B-33). 


3.  Implications 


In  this  report,  we  derive  equations  describing  the  scattered  fields  from  a 
spherical  scatterer  located  eccentrically  within  a  host  sphere.  The  source  is 
also  arbitrarily  located  within  the  host  sphere,  and  can  be  chosen  to  simulate 
a  P  wave,  S  wave,  or  any  linear  combination  of  the  two.  We  can  use  this  theory 
in  seismology  to  predict  the  scattering  interaction  between  a  scattering  body 
and  the  Earth’s  surface. 

The  scattered  fields  are  complete,  and  exact  solutions  for  the  elastic  wave- 
field  produced  by  the  spherical  scatterer  are  expressed.  Because  the  solution 
presented  here  is  almost  entirely  analytical,  it  avoids  the  need  for  much  of 
the  computationally  intensive  algorithms,  such  as  finite-difference  schemes, 
that  would  normally  be  required  to  solve  this  problem.  Also,  because  we 
have  made  very  few  assumptions  concerning  the  boundary  conditions,  the 
analytical  solution  given  here  should  be  well  behaved. 

This  theory  provides  a  basis  that  should  prove  useful  in  understanding  and 
predicting  seismic  scattering.  Using  techniques  developed  in  electromagnetic 
scattering  (Videen  et  al,  1995),  it  should  be  possible  to  extend  this  theory 
to  describe  the  scattering  from  an  arbitrary  body  rather  than  a  spherical 
one.  Hence,  this  theory  is  an  important  first  step  in  demonstrating  a  method 
for  analytically  deriving  scattered  elastic  wave  fields  from  anomalies  such  as 
magma  chambers  or  subducting  slabs  near  the  Earth’s  surface. 
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Appendix  A.  Boundary  Coefficients 


We  can  find  the  coefficients  that  satisfy  the  boundary  conditions  at  the  scat¬ 
tering  sphere  outer  surface  (r2  =  a')  by  first  eliminating  the  dependence  on 
the  internal  field  coefficients.  Performing  some  algebraic  manipulation  and 
taking  advantage  of  the  differential  equation  defining  spherical  Bessel  func¬ 
tions, 

(r)  +  2^4  {r)  +  -  n  (n  1)]  Zn  (r)  =  0  ,  (A-1) 


we  can  put  the  coefficients  that  satisfy  the  boundary  conditions  on  the  scat¬ 
tering  sphere  into  the  following  form: 


^  nm 

A\Bl  -  A\B\ 
B\Bl-B\B\  ’ 

(A-2) 

_ 

^nm 

A\Bl  -  AlB\ 

B\Bl-  B\Bl  ’ 

(A-3) 

^nm 

4(1,1) 

4(3,1)  ’ 

(A-4) 

0®^  = 
Vnm 

A\Bl  -  A\B\ 

B\Bl-B\Bl  ’ 

(A-5) 

II 

A\Bl  -  AlB\ 

B\Bl-B\Bl  ’ 

(A-6) 

oil  — 
^nm 

C\Bl  -  ClBl 
BlBi-BiB!  ’ 

(A-7) 

pl3  _ 

CfBl-CiBl 

B\Bl-  B\B\  ’ 

(A-8) 

p22  __ 

■^nm 

4(4,1) 

4(3,1)  ’ 

(A-9) 

p31  _ 

^nm 

C\Bl  -  ClB\ 

B\Bl  -  B\Bl  ’ 

(A-10) 

o33  _ 

^nm 

ClB\-ClB\ 

B\Bl  -  B\B\  ’ 

(A-11) 

where 


and 


=  {A';Dl-A^Dl){DiDl-DiDl) 

-  (A^Dl  -  AlDtj  (Dil^  -  D‘d!) 

(A-12) 

An 

^2 

=  (A^DI  -  A^D!)  (D‘Di  -  DlDl) 

-  [A^Dl  -  AlDf)  {d\i4  -  D’DJ) 

(A-13) 

B'l 

=  (B^DI  -  B^Df)  (D‘D^  -  D}Dl) 

-  (BJDJ  -  b;dI)  (dIdI  -  DiD\] 

1  {A-14) 

=  {b^dI  -  B^Df')  (d‘D|  -  D^Dl) 

-  (BJB|  -  BJD’)  (d;dI  -  B^Df; 

1  (A-15) 

=  (CfBl  -  CJDf)  (D^D!  -  D‘Di) 

-  {CSDl  -  c;dI)  (bJdI  -  B‘D?) 

\  (A-16) 

nn 

=  (q-Dl  -  C^Df)  (D>B’  -  D‘Dl) 

-  (C^Dl  -  CJB|)  {DlDl  -  DlLkl] 

1  (A-17) 

=  i/2j— 1  (1?  ^a?  ^5  ^  )  j 

(A-18) 

^  H2j  X,  a  ) , 

{A-19) 

B] 

==  Il2j—l  (3,  ko,^  X^  ix^ci  ^  ^ 

(A-20) 

B] 

=  H2j  (3,  A,  //,  a') , 

(A-21) 

^2j—l  (4?  ^a5  A,  (X^  ^  1 

(A-22) 

”  (4j  k^^  A,  /i»  (2  ) , 

(A-23) 

D] 

=  -f^2j-l  (I7  ^a’  ^0  ’ 

(A-24) 

=  (1,  A', /z',  a')  , 

(A-25) 

Jniit  j')  —  Hfi (z ,  kp ,  X,  fJ,,  (I  ) 


r(») 

Zn 
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Hi{p,k,X,n,a)  =  2p 


(A-26) 


H2{p,k,X,fi,a) 

H3{p,k,X,p,  a) 
H4{p,k,X,p,a) 


=  2ii(n*  +  n-  _  2itzlF'>'(ka)  -  likaz^S'’ (ka)  ,  (A-27) 

ka 


2pn{n  +  1) 


ka 


—  ka{X  +  2p) 


z^\ka)  -  4pzl^^'{ka)  , 


=  pn{n  +  l)Hi{p,k,a) 


(A-28) 


The  coefficients  that  satisfy  the  boundary  conditions  on  the  host-sphere  outer 
surface  (ri  =  a)  are 


?ii 


?13 


?22 


ySl 


:-33 


kg^  A,  p,  q) ^^2(4,  fcjg)  A)  p,  a)  H\{3.  kg,  A,  p,  n) ^4(4,  fcg,  A,  p,  q)  29) 

i/i(4,  A:a,  A,  p,  a)H4{4,  kp,  X,  p,  a)  -  ^3(4,  ka,  A,  p,  a)H2{4,  kp,  A,  p,  a)  ^ 

114(8,  kp,  X,  p,  a')  112(4,  kp,  X,  p,  fl)  112(8,  kp,  X,  p,  a')H4(4,  kp,  X,  p,  o)  30) 

Hi(4,  ka,  X,  p,  a)H4(4,  kp,  X,  p,  a)  -  Hz(4,  ka,  X,  p,  a)H2(4,  kp,  X,  p,  a)  ^ 


Hi(8,kp,X,p,a) 

7^1  (4,  kp,  X,  p,  a) 

llz(^,  ^cc,  Ai  /^5  a'^H 1(4,  ka,  X,  p,  a^  Hi (3^  ka,  X,  p,  a^H^ (4,  ka,  X,  p,  a) 
113(4,  ka,  X,  p,  a)H2(4,  kp,  X,  p,  a)  -  Hi(4,  ka,  X,  p,  a)H4(4,  kp,  X,  p,  a) 
114(8,  kp,  X,  p,  a'^Hi (4,  ka,  X,  p,  a)  112(8,  kp,  X,  p,  a) (4,  ka,  X,  p,  a) 
ff3(4,  ka,  X,  p,  a)H2(4,  kp,  X,  p,  a)  -  Hi(4,  ka,  X,  p,  a)H4(4,  kp,  X,  p,  a) 


(A-31) 

(A-32) 

(A-33) 
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Appendix  B.  Vector  Translations 


Stein  (1961)  and  Cruzan  (1962)  derived  translation  addition  theorems  for 
vector  spherical  wave  functions  that  can  be  used  to  express  the  displacement 
fields  in  one  coordinate  system  in  terms  of  displacement  fields  in  another 
coordinate  system.  These  results  are  also  given  by  Ben-Menahem  and  Singh 
(1981).  We  use  their  results  to  find  relationships  between  the  coefficients 
and  and  and  and  For  a  translation  along  the  z-axis 

with  no  rotation,  the  vector  spherical  harmonics  are  related  by 


oo 

T  (p)  _  {?) 

^nm,2  /  V  ’ 

(B-l) 

n'=0 

CX) 

(B-2) 

n'=0 

IVjCp)  _  ^  R(’^’0Tvr(9) 

(B-3) 

n'=0 


If  p  =  3  or  4,  then  z  =  1  in  the  region  where  ri  <  |d| ,  and  ^  =  p  in  the  region 
where  n  >  \d\ .  Recurrence  relations  for  the  scalar  translation  coefficients 
were  derived  by  Bobbert  and  Vlieger  (1986),  which  greatly  facilitate 
calculations.  The  results  are 


=  vwnzi^){kd), 

(B-4) 

c[y„,  =  -V2n'  +  lzll){kd)  , 

(B-5) 

2n'+l 


2iTi  +  1 


(B-6) 


v/(r!  -  m  +  l)(n  +  m)(2n'  +  -m+l){n'  +  m)(2n'  +  -  (B-7) 

kd. 


kd 


{n'  -m  +  2)(n'  -  m  +  1) 

^n,n'+l 


(2n'  +  3) 


(n'  +  m){n'  +  m-  1) 
(2n'  -  1) 
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From  the  scalar  translation  coelRcients,  the  vector  translation  coefficients  can 
be  easily  derived: 


kd 


'n,n' 


n 


{n'  -m  +  l)(n'  +  m  +  1) 


'  +  1^  (2n'  +  l)(2n'  +  3) 


n,n'H-l 


(B-8) 


kd  (n'  — m)(n'  +  m)  (jni) 
'^\  (2n'  +  l)(2n'-l)  ’ 


-imkd 


a 


(m,i) 


n'(n'  +  l) 


(B-9) 


Prom  these  equations,  we  see  that 


A{m,p,q)  _  A-m,p,q) 
^n,n'  ^n,n'  ’ 


(B-IO) 


r>(»”.P.9)  _  r>{-m,p,q) 


(B-U) 


^n,n'  ^n,n' 


(B-12) 


Since  uj„  =  u^,  =  u^,  and  u‘„,  =  uL,  we  can  derive  expressions 

relating  the  displacement  coefficients  using  equations  (B-1)  to  (B-3); 


aill  = 

oo 

n'=:0 

(B-13) 

0!^nl  = 

oo 

n'=0 

(B-14) 

a®  = 

nm 

oo 

E  ’(-M) 

n'=0 

(B-15) 

fW  - 
Jnm 

oo 

n'=0 

(B-16) 
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00 

E  (M)  +  h‘ilB'S‘^i\k,d) . 

n'=0 

(B-17) 

/•(3)  _ 

Jnm 

oo 

n'=0 

(B-18) 

= 

onm 

oo 

n'=0 
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(B-20) 


(7(2) 

^nm 

n'=0 

cfnm 

oo 

—  V  r(^^ 

n'=0 

^(2)  r>(jT^,l)/ 


(B-21) 


Finally,  we  note  that  the  scalar  translation  coefficients  used  in  deriving 
the  vector  translation  coefficients  and  are,  in  general,  different 

from  the  vector  translation  coefficients  used  in  the  translation  of  the 

P  waves  (La).  since  in  general  the  P  and  S  waves  have  different  spatial 
frequencies.  In  addition,  the  translations  for  the  source  coefficients  and 
are  in  opposite  directions. 

The  translation  from  the  source  coordinate  system  to  the  rotated  host  coor¬ 
dinate  system  is  similar: 


i(i)  = 

^nm 

oo 

n'=0 

(B-22) 

t(2)  = 

oo 

n'=0 

(B-23) 

*(3)  = 

oo 

(B-24) 

n'=0 


Using  the  results  of  Stein  (1960),  we  present  the  coefficients  in  a  rotated 
coordinate  system.  The  vector  harmonics  can  be  expanded  as 


(p) 

'nmj 

—  n  ("’’”)  T  (p) 

2^  ^nm\i  ’ 

(B-25) 

m'=0 

(p) 

nmj 

oo 

(B-26) 

m'=0 

■(p) 

nmj 

oo 

—  n(«.»p)Aj{9) 
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where 


21 


(B-28) 


exp  [i  {m'a  +  m'y)] 
n  +  m  '' 


{n  +  m')\  (n  —  m')\ 
(n  +  m)\  {n  —  m)!  J 


11/2 
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E 

(7 

[sin(/?/2)] 


n  —  m—aj\ 

2n—2a—m'—m 


(-1) 


n+m-cr 


cos  (/?/2)] 


2(7+m'+m 


(B-29) 


and  a,  (3,  and  7  are  Euler  angles  using  the  convention  of  Edmonds  (1957). 
The  coefficients  in  the  rotated  coordinate  system,  CLnm,ii  can  be  expressed  in 
the  unrotated  coordinate  system,  anmj^  as 


a 


nmj 


=  Eo 

m'=0 


(n,m) 


(B-30) 


Therefore,  the  source  coefficients  in  the  host  coordinate  system  for  the  point 
source  can  be  expressed  as 


gd)  = 


e(2)  = 


gO)  = 


00  00 

n'=0  n'=0 
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m'=0 
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